Gravitational baryogenesis in Gauss-Bonnet braneworld cosmology 
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The mechanism of gravitational baryogenesis, based on the CPT-violating gravitational inter- 
action between the derivative of the Ricci scalar curvature and the baryon- number current, is in- 
vestigated in the context of the Gauss-Bonnet braneworld cosmology. We study the constraints 
on the fundamental five-dimensional gravity scale, the effective scale of B- violation and the decou- 
pling temperature, for the above mechanism to generate an acceptable baryon asymmetry during 
the radiation-dominated era. The scenario of gravitational leptogenesis, where the lepton-number 
violating interactions are associated with the neutrino mass seesaw operator, is also considered. 
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INTRODUCTION 



The origin of the baryon asymmetry is an outstand- 
ing problem in particle physics and cosmology. Suf- 
ficient conditions for baryogenesis are the violation of 
baryon number, the violation of C and CP symmetries 
and the existence of nonequilibrium processes pj. Al- 
ternatively, if CPT and baryon number are violated, a 
baryon asymmetry could arise even in thermal equilib- 
rium 0, Q. In Ref. 0, the effects on baryogenesis of 
certain CPT-violating terms arising in a string-based 
framework were investigated and it was shown that a 
large baryon asymmetry could be produced at the grand- 
unified scale. Recently, a new baryogenesis mechanism, 
where the baryon asymmetry is generated via a dynami- 
cal breaking of CPT while maintaining thermal equilib- 
rium, was proposed in Ref. [j] . The crucial ingredient is 
an interaction between the derivative of the Ricci scalar 
curvature TZ and the baryon-number (B) current J M (or 
any current that leads to a net B — L charge in equilib- 
rium, where L is the lepton number, so that the asymme- 
try will not be erased by the electroweak anomaly 6J): 



1 

W 



d 4 arv / =0(d M 7e)J'\ 



(1) 



where M» characterizes the scale of the interaction in the 
effective theory. Such an operator is expected to arise in 
the low-energy effective field theory of quantum gravity 
or in supergravity theories from a higher-dimensional op- 
erator 0. 

The interaction in Eq. violates CP and, in an ex- 
panding universe, it also dynamically breaks CPT. If one 
requires the existence of P-violating processes in thermal 
equilibrium, then a net baryon asymmetry can survive 



after their decoupling at a temperature Td @ , 



ns TZ 



(2) 



For this mechanism to work, a nonvanishing time deriva- 
tive TZ 7^ is necessary. 1 Although in an expanding 
universe TZ ~ H 2 is nonzero in four-dimensional gen- 
eral relativity (GR), its time derivative TZ = during 
the radiation-dominated (RD) epoch. It turns out, how- 
ever, that TZ 7^ can be easily realized in the braneworld 
scenario, which suggests that higher-dimensional gravity 
effects can offer a novel way to generate a baryon asym- 
metry through the dynamics of spacetime [t| . 

In the past few years, stimulated by the development of 
string theory, the braneworld id eas, and particularly, the 
Randall-Sundrum (RS) model [T(|, have been actively 
investigated. The RS braneworld cosmology is based on 
the five-dimensional Einstein-Hilbert action; at high en- 
ergies, it is expected that this action will acquire quan- 
tum corrections, in the form of higher-order curvature 
invariants in the bulk action. String theory and holog- 
raphy indicate that such terms arise in the action at the 
perturbative level. In five dimensions, the Gauss-Bonnet 
(GB) invariant has special properties: it represents the 
unique combination that leads to second-order gravita- 
tional field equations linear in the second derivatives and 
is ghost- free [Til, u2, LL21 03 • Moreover, the graviton zero 
mode remains localized in the GB braneworld |l5j and 
deviations from Newton's law at low energies are less 
pronounced than in the RS case ^(J . 

In this paper we examine gravitational baryogenesis 
in the context of GB braneworld cosmology. The case 
when the GB contribution is absent and cosmology is of 
RS type is also considered. We show that the observed 
baryon-to-cntropy ratio can be successfully explained in 
both frameworks. The possibility that B — L violation 
is associated with the neutrino mass seesaw operator is 
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1 For a more general form of the derivative coupling of the Ricci 
scalar to ordinary matter, C oc d/j,f(lZ) J M , see Ref. [jj. 
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also studied. In the latter case, the limits coming from 
low-energy neutrino physics when combined with the GB 
inflationary constraints allow us to put bounds on the 
fundamental scale of gravity the effective scale of B- 
violation and the decoupling temperature, which are re- 
quired to generate an acceptable baryon asymmetry in 
the gravitational leptogenesis scenario. 



II. GAUSS-BONNET BRANEWORLD 

The five-dimensional bulk action for the Gauss-Bonnet 
braneworld scenario is given by 

S = W-2 i d 5 x^-®g [-2A 5 + R 
Zk 5 J 

+ a(R 2 - AR ab R ab + R abcd R abcd )] 



d A x \f^g A + S n 



(3) 



branc 



where a > is the GB coupling, A > is the brane 
tension, A5 < is the bulk cosmological constant and 
S'mat denotes the matter action. The fundamental energy 
scale of gravity is the 5D scale M5 with n\ = 8%/M§, and 
M 4 is the 4D Planck scale with k 2 = 8n/M 2 . 

The GB term may be viewed as the lowest-order 
stringy correction to the 5D Einstein-Hilbert action with 
a <C £ 2 , where t is the bulk curvature scale, \R\ ~ l/l 2 . 
The Randall-Sundrum type models are recovered for 
a = 0. Moreover, for an anti-de Sitter bulk, it follows 
that A 5 = —3p 2 (2 — /?), where p = l/l is the extra- 
dimensional energy scale and 



(3 = 4ap 2 <C 1 . 



(4) 



Imposing a Zi symmetry across the brane in an anti- 
de Sitter bulk and assuming that a perfect fluid matter 
source is confined to the brane, one gets the modified 
Friedmann equation |12L Il7| 



4(p + \) = 2/iJl + 



IP 

J 2 



3-/3 + 2/3- 



P 2 



This can be rewritten in the useful form 1' 







(l-/3)cosh(^) -1 



(5) 



(6) 



where x is a dimensionless measure of the energy density 
p on the brane such that 



with 



p + A = m a sinh \ , 



8p 2 (l - pf 



P4 



(7) 



(8) 



the characteristic GB energy scale. The GB high-energy 
regime (sinhx ^> 1) corresponds then to p 3> m^- Notice 



also that we must have m a > m\ = A 1 / 4 , where m\ is 
the characteristic RS energy scale. This in turn implies 

/3< 0.15 m. 

The requirement that one should recover general rela- 
tivity at low energies leads to the relation [l2l IT^ | 



/' 



(9) 



Since /3 <C 1, we have p, w Af|/M|. Furthermore, the 
brane tension is fine-tuned to achieve a zero cosmological 
constant on the brane 1191: 



K 2 5 X 



2p(3-(3) 



(10) 



Expanding Eq. J^J in %, we find three reg imes for the 
dynamical history of the brane universe |18j 



p > m a 



to 4 » p » m\ 



p <C m x 



ll 2 K 2 
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4/3 
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(GB), 

(RS), 
(GR). 



(11) 
(12) 
(13) 



Eqs. (|ll|) - (|13|l are considerably simpler than the full 
Friedmann equation and in many practical cases one of 
the three regimes can be assumed. In this case, it is useful 
to consider a single patch with the effective Friedmann 
equation ppf 



H 2 



(14) 



where q = 1, 2, 2/3 for GR, RS and GB regimes, re- 
spectively. For each regime, the coefficients f3 q > are 
determined in accordance with Eqs. (|11[1 - H13|) . 



III. GRAVITATIONAL BARYOGENESIS 

Let us now consider the gravitational baryogenesis 
mechanism in the GB braneworld. In an expanding uni- 
verse, the interaction term in Eq. gives rise to an 
effective chemical potential p b ~ 1Z/M 2 for baryons. In 
thermal equilibrium, the net baryon-number density does 
not vanish as long as pb ^ 0, and for mb, pb *C T one 
has H3 



n B 



9b 
6 



PbT 2 



(15) 



where gb is the number of intrinsic degrees of freedom 
of the baryon. During the RD epoch, this leads to a 
baryon-to-entropy ratio given by 



Td 



with 



15ffb 



(16) 



(17) 



3 



We have used s = 27r 2 <7* s T 3 /45 for the entropy density, 
where g* s is the total number of degrees of freedom which 
contribute to the entropy of the universe. 

The Ricci scalar in the Friedman-Robertson- Walker 
brane is defined in terms of the expansion law as 

K = -6(H + 2H 2 ). (18) 

Using the full GB Fricdmann Eq. © , we get for the time 
derivative of the Ricci scalar, 



n 

ri(x) 



4fx 2 {l + w)(l-/3) Hp 



6 + ; 
3(1 + w' 



(3 mi 



w) 



sinh 



coshx 
2x 



[n(x) +r 2 (x)} , 



mi coshx 



co.dl ( 



— sinh | — I tanh v 
2 V 3 ' A 



(19) 



where w = p/p is the equation of state. During the RD 
era, w = 1/3, ri(x) = and Eq. 1)19(1 simplifies to 

7t = _32 M 2 (l~ / 3) iJp 2 



(3 m| cosh 2 x 

- cosh | — | — sinh ( — I tanh y 
3 V 3 / V 3 ' 



(20) 



Combining Eqs. IjRjfl and H2U|) . it is then possible to 
compute the decoupling temperature, Tjj, required to 
produce an acceptable baryon-to-entropy ratio ns/s. In 
Fig. we plot this temperature as a function of the scale 
Af* for different values of the GB coupling (3 and the fun- 
damental gravity scale M5, assuming the observed value 
n B /s ~ 9 x 10 -11 [23. We have used the fact that in the 
RD era the energy density is p = 7r 2 g >t T 4 /30, where g* 
is the total number of relativistic degrees of freedom. In 
the standard model ~ g* s ~ 100 above the electroweak 
scale and, assuming g^ ~ 0(1), one gets c ~ 0(1O~ 2 ). 

In order to establish the transition between the differ- 
ent regimes, we can consider the simplified expansion law 
given by Eq. I|14|). In this case, 



K = -48q(q~l)H 3 



(21) 



and the baryon asymmetry reads as 



2 x3 9 /2 6,-1 



If the decoupling occurs in the RS regime, where q = 2 
and f3 q = (K4/6A) 1 / 2 , we obtain 



2 1**2 Ji/r9\ i/ 11 



To ^ 3.2 x 10~ 2 (M 2 M|) 



(23) 



for (3 -C 1. For the decoupling to occur in this regime, it 
is required that p(Tjj) > m\ , which implies in turn that 
M* > M^ s with 



/MX 11 / 4 
/ !' 10 5 (^] M 5 



(24) 



On the other hand, in the GB regime, where q = 2/3 and 

1 /3 

f3 q = (/i 2 K§/4/3) , the decoupling temperature is given 
by 



Tn ~ 1.6 X 10 



_ 4 (/3M|M ( 
M 5 



4x1/3 
4 J 



(25) 



so that the transition from the GB to the RS regime, 
defined by the condition p(Tb) = m 4 , occurs for 



GB 



5.9 x 10 4 /Af 5 \ n/4 



A I: 



5 • 



(26) 



The transition values M^ 5 and are represented in 

Fig. G] by the two vertical dashed lines. We see that the 
decoupling of the baryon-number violating interactions 
can generally occur in either (GB, RS or GR) regime. 
We also notice that the RS transition region shrinks as 
the Gauss-Bonnet coupling (3 increases. 

An important constraint on the decoupling tempera- 
ture comes from reheating and inflation. Obviously, we 
must require Tjj < T r h < Mj, where T r h is the reheating 
temperature (at which the universe becomes radiation 
dominated) and Mj is the inflation scale. In the con- 
ventional scenario, reheating occurs as the inflaton field 
oscillates around its minimum and decays into matter. 
In this case, T r h crucially depends on the details of the 
inflaton coupling to matter. Here we restrict our discus- 
sion to the inflation scale Mj. Assuming that inflation 
occurs in the high-energy (RS or GB) regime of the the- 
ory and that it is driven by a quadratic potential, we 
find Mi w m a sinh 1/,4 x e , where \& is evaluated at the 
end of inflation (see the Appendix for details) . In Fig. ^ 
we have plotted Afj for different values of (3 and A/5. We 
note that for Af* < M5 the constraint Tr> < Mi is always 
verified. 

We have also performed a random scan of the parame- 
ter space (Af*, Mg,/3) in order to find the allowed region 
for the gravitational baryogenesis mechanism considered 
here to generate an observationally acceptable ub/s with 
T D < Mj < M 5 . The results are presented in Fig. A 
similar analysis was done for the case when the GB terms 
are absent, i.e. (3 = 0, and braneworld cosmology is of 
RS type (see Fig. • We notice that in the RS case the 
gravity scale M5 can take considerably lower values (cf. 
Fig. (3a)), only constrained to be larger than 10 5 TeV, 
if one requires the theory to reduce to Newtonian grav- 
ity on scales larger than 1 mm. The above bound yields 
T D > 10 5 GeV. 

Up to now we have not taken into account pos- 
sible effects which could dilute the baryon asymme- 
try generated by the mechanism described above. It 
is well known that electroweak sphaleron transitions, 
which are unsuppressed at temperatures above the elec- 
troweak phase transition, are a potential source of dilu- 
tion 0| . Sphaleron- induced baryon-asymmetry dilution 
occurs when B — L vanishes. In this case, the B and 
L-number densities will be typically diluted by a factor 
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FIG. 1: The decoupling temperature To and the scale Mb of B — L violation as functions of M* for different values of the 
Gauss-Bonnet coupling f3 and the fundamental scale of gravity Ms. The vertical dashed lines delimit the transitions between the 
GB, RS and GR regimes. The inflation scale Mi (dashed line) and the decoupling temperatures (dot-dashed lines) associated 
with the neutrino mass seesaw operator in the case of quasidegenerate (T® D ) and hierarchical (T^ 11 ) neutrino masses are also 
shown. 



0.02 m 2 /T 2 ph Q, where m T is the t lepton mass and T sph 
is the sphaleron freeze-out temperature. Assuming T sp h 
to be the electroweak scale, one finds that the baryon 
asymmetry is diluted by a factor of about 10~ 6 . Hence, 
according to Eq. (J2J), the scale M* would have to be, in 
this case, smaller by a factor of 10~ 3 to reproduce the 
correct value of ns/s via the gravitational interaction of 
Eq. QJ. On the other hand, if B — L ^ 0, essentially no 
sphaleron dilution occurs. In the latter case, the baryon 
asymmetry generated will remain after the decoupling of 
the (B — L)-violating interactions. An example of this 
possibility will be presented in the next section. 



IV. GRAVITATIONAL LEPTOGENESIS 

In the standard model of electroweak interactions, the 
B — L symmetry is exactly conserved. This symmetry is 
however violated in many of its extensions. In general, it 
is possible that the B- violating interactions are generated 



by an operator Ob of mass dimension D = 4 + n. The 
rate of such interactions is Tb ~ T 2n+1 /M B n , where Mb 
is the mass scale associated with the operator Ob- In 
the standard electroweak model the lowest-dimensional 
operator that violates B—L is the dimension five operator 

£f = jgt£4><f> + B.c, (27) 

where I and cf> are the left-handed lepton and Higgs dou- 
blets, respectively; M is the scale of new physics which 
induces B — L violation. This interaction represents a 
typical term that gives rise to the seesaw mechanism and 
is responsible for the light neutrino masses ~ v 2 /M, 
v ~ 174 GeV. In the early universe the L-violating rate 
induced by the interaction l|27|) is [2^] 

T 3 10 v 2 

r ^ = M|' (E-f) 1/2 ' (28) 

The decoupling of the (B — L)-violating processes oc- 
curs when falls below the Hubble rate, i.e. when 
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FIG. 2: The parameter space that generates an observationally acceptable baryon asymmetry during the radiation era in 
Gauss-Bonnet braneworld cosmology with To < Mi < Ms. We take ub/s — 9 x 10 -11 and assume that inflation is driven 
by a quadratic potential in the high-energy (RS or GB) regime of the theory. The black-dotted region is obtained by taking 
into account the constraint Mb < Mi for the mass scale associated with a dimension five B-violating operator. The horizontal 
dot-dashed lines correspond to the GR decoupling temperatures T® D and T^ 1 (panel (c)) and the scales Mg D and M§ ! (panel 
(d)) associated with the neutrino mass seesaw operator. 



Mb — [Tp/if (Td)] 1 / 2 . In the gravitational baryogene- 
sis scenario considered here, the decoupling temperature 
that produces an acceptable baryon asymmetry is fixed 
by Eq. 1|16[) , and it determines the required scale of B — L 
violation. In Fig. ^ we have plotted Mb as a function of 
M* for different values of f3 and M5. We notice that the 
requirement To < Mb < M5 imposes an upper bound 
on M*. We find M* < 10 16 GeV. 



If the scale Mb is associated with the neutrino mass 
seesaw operator, as in Eq. I|28[) . then the value of this 
scale will be fixed by the light neutrino mass spec- 
trum. The current cosmological limit coming from the 
Wilkinson Microwave Anisotropy Probe (WMAP) im- 
plies ^2 nii < 0.69 eV p^ . If neutrinos are quasidegen- 
erate (QD) in mass, the above limit requires m\ ~ mi ~ 



to 3 ~ 0.23 eV. In this case, 

Mg D w 7.6 x 10 14 GeV . (29) 

Instead, if neutrinos masses are hierarchical (HI) with 
mi ~ <C to 2 "C 1TI3, then to 2 ~ (Am^) 1 ' 2 and 
m.3 ~ (Am 2 ^) 1 / 2 , where the squared mass differences 
measured in solar and atmospheric neutrino oscillation 
experiments are Ato 2 o1 ~ 8.1 x 10~ 5 eV 2 and Am 2 tm ~ 
2.2 x 10~ 3 eV 2 0, respectively. In the latter case, 

Mg 1 s» 6.3 x 10 15 GeV. (30) 

Eqs. I|29(l and (|30(l yield a decoupling temperature in the 
range 

T® D < T D < T^ 1 . (31) 
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FIG. 3: The parameter space that generates the observed baryon-to-entropy ratio tib/s during the radiation era of Randall- 
Sundrum braneworld cosmology (j3 = 0). In this case, Mi m 5 x l0~ 2 Ms. For other details, see caption of Fig. [5] 



Clearly, the specific values of T® D and T^ 1 depend on 
whether the decoupling of B — L violation occurs in GB, 
RS or GR regime and, thus, on the values of the Gauss- 
Bonnet coupling (3 and the fundamental scale M$. In 
standard cosmology with H(Te>) ~ 1.66 gl^T^/M^ , one 
finds T® D « 7.9 x 10 11 GeV and T? 1 rs 5.5 x 10 13 GeV. 
Some other examples are presented in Fig. ^ While 
in Fig. (la) and (lb), the decoupling corresponding to 
T® D and T^ 1 (horizontal dot-dashed lines) occurs in 
standard cosmology, in Fig. (lc) and (Id) such decou- 
pling takes place in the high-energy Gauss-Bonnet regime 
for the case of hierarchical neutrinos. We also remark 
that for gravitational leptogenesis to be successful we 
must require To > T^ D , which implies the lower bound 
M* > 100 GeV for f3 < 0.1, as can be seen from the 
figure. 

Let us now consider the inflation bound. Since To < 
Mb, the requirement Mb < Mi is more stringent in this 
case. We find that this bound strongly constrains the 
scale of B — L violation and, consequently, the mecha- 
nism of gravitational leptogenesis. For instance, it can be 



seen that for the case presented in Fig. (Id), the above 
constraint implies the bound Tjj < T® D and, therefore, 
the leptogenesis mechanism cannot generate the required 
baryon asymmetry. The allowed region for gravitational 
leptogenesis is presented in Fig. [21 (black dots). The hori- 
zontal dot-dashed lines correspond to the GR decoupling 
temperatures T® D and T^ 1 (Fig. (2c)) and the scales 
M% D and M§ 1 (Fig. (2d)) associated with the neutrino 
mass seesaw operator. We conclude that 

10 15 GeV< M 5 < 10 17 GeV, 

10 2 GeV < M* < 10 10 GeV . (32) 

One can compare the above results with the ones that 
are obtained in the case when braneworld cosmology is 
of RS type, i.e. when (3 — 0. The allowed range of values 
for the parameters is shown in Fig. [21 We notice that 
a successful gravitational leptogenesis in RS cosmology 
requires 

10 16 GeV < M 5 < 10 17 GeV , 

10 2 GeV < M* < 10 6 GeV . (33) 
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CONCLUSION 



that inflation is driven by the simple quadratic potential 



In this work we have considered the possibility that 
the observed baryon asymmetry arises via the spacetime 
dynamics of Gauss-Bonnet braneworld cosmology. The 
framework presented here is based on the CPT- violating 
gravitational interaction between the derivative of the 
Ricci scalar curvature and the B (or B — L) current 
We have shown that it is possible to generate the correct 
magnitude of the baryon asymmetry in different cosmo- 
logical scenarios, depending on whether the decoupling of 
the B- or (B — i)-violating interactions occurs in stan- 
dard cosmology or in the high-energy Randall-Sundrum 
or Gauss-Bonnet braneworld regimes. 

We have also studied the case when baryogenesis oc- 
curs via leptogenesis, and the B — L current is associ- 
ated with the neutrino mass seesaw operator. In this 
framework, the produced ub-l asymmetry will be con- 
verted to a baryon asymmetry, once sphalcron transi- 
tions enter thermal equilibrium. We have seen that 
for this scenario to be viable, a rather high fundamen- 
tal scale of gravity M5 is required (cf. Eq. iJjUJl), as 
well as an effective interaction scale M* above the elec- 
troweak scale but below 10 10 GeV. At this point it is 
worth noticing that, although in four-dimensional grav- 
ity it is natural to expect M* of the order of the Planck 
mass M4, this may not be necessarily the case. For 
instance, M* ~ (Mr A/4) 1 / 2 could be possible, if the 
right-handed neutrino Majorana mass Mr softly violates 
baryon number ,5]. Moreover, this scale can be much 
lower, if the effective four-dimensional theory comes from 
a higher-dimensional theory. Indeed, AdS/CFT corre- 
spondence [2^| and braneworld holography 0, [2(| in- 
dicate that interaction terms such as given by Eq. (JTJ 
are expected in the effective action on the brane with 
M* ~ 1/t ~ Ml/ Ml. If this is the case, the bound M 5 ~ 
10 15 - 10 17 GeV would then imply M* ~ 10 7 - 10 13 GeV, 
well within the range allowed by gravitational baryogen- 
esis and leptogenesis in the high-energy GB regime. 
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APPENDIX: INFLATION IN A GAUSS-BONNET 
BRANEWORLD 

In this Appendix we briefly review inflation in GB 
brane cosmology 0, 0, |53 • For simplicity we assume 



V(<j)) = V 4> 2 - (A.l) 
We are interested in slow-roll inflation which occurs 
in the GB or RS high-energy regime. In this case, V » 
p ^> A and Eq. J7J implies V « sinhx- Moreover, the 
bound V > m\ together with the quantum gravity upper 
limit V < M| imply 



— -j < sinhx < ( — - 
m a \m a 



The slow-roll parameters e and 77 are given by 



(A.2) 



16AVb 

27K?m8 



fix) > V 



8XVn 1 



where 



9K|m8 g(x) 



2 X 



(A.3) 



fix) =9 2 (x)tanhxsinh — , (A.4) 



,'/(\! : (^1 - 1 



The number of e-folds of inflation is given by 



v, =^tt f X gix)coth x d x = ^f j (x) 



(A.5) 



where 



Iix) = 9ix) ~ ^ 



8PV 



1 + 3<7(X) 



(A.6) 



(A.7) 



X* is evaluated when cosmological scales leave the hori- 
zon and Xe is evaluated at the end of inflation, when 
max{e, 77} = 1. 

The amplitude of scalar perturbations is 



Ai 



35/2 



ff 3 (x*) 



l&^VoX 1 / 2 ^/ 2 sinhx* 



(A. 



Using the COBE normalized value As — 2 x 10 5 for the 
density perturbations and 55 < N+ < 65, we can obtain 
the scale of inflation Mi = V 1 ^ 4 (4> e ) w m a sinh 1 / 4 x e ■ 
This scale is plotted in Fig^for given values of (3 and M5, 
taking JV* = 60. A more complete analysis is presented 
in Fig. (2a) and (2b). Notice that, for consistency, one 
should require m\ < Mj < M5 . 
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